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The linear potential binding a quark and an antiquark in mesons is generalized to baryons and multiquark 
configurations as the minimal length of flux tubes neutralizing the color, in units of the string tension. For 
tetraquark systems, i.e., two quarks and two antiquarks, this involves the two possible quark-antiquark pairings, 
and the Steiner tree linking the quarks to the antiquarks. A novel inequality for this potential demonstrates 
rigorously that within this model the tetraquark is stable in the limit of large quark-to-antiquark mass ratio. 
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The quark-antiquark confinement in ordinary mesons is often described by a linear potential V2 = r, in units where the string 
tension is set to unity. For a given interquark separation r, it can be interpreted as the minimal gluon energy if the field is 
localized in a flux tube of constant section linking the quark to the antiquark. 

The natural extension to describe the confinement of three quarks in a baryon is the so-called F-shape potential 

V3(wi,i'2,W3) = niin(di +^2 + ^3) , (1) 

where di is the distance of the quark located at Vi {i = 1, 2, 3) to a junction s whose location is adjusted to minimize V3. 
This potential has been proposed in Refs. ||Tl|2][3]|4l|5]|6l|2l, among others. It has been used, e.g., in Refs. |[8]|9l for studying 
the spectroscopy of baryons. See, also (TQ\. The optimization in (p^ corresponds to the well-known problem of Fermat and 
Torricelli to link three points with the minimal network. See Fig. |l(a)| 

We now turn to the tetraquark systems {Q, Q, q, q), with the notation {vi, V2,V3, V4) for the locations, and (Af, M, m, m) for 
the masses which will be used shortly. The potential is assumed to be (with dij = \\viVj\\) 



U = min{(ii3 -|- (^24, du + ^23, Vi] ■ 

Vi = min ( ll-wisill + ||w2Si|| + ||siS2|| + IIS2W3II + ||s2W4|| ) • 

Si,S2 



The first two terms of U describe the two possible quark-antiquark links, and their minimum is sometimes referred to as the 
"flip-flop" model, schematically pictured in Fig. 1(b) It was introduced by Lenz et al. ifTTl . who used, however, a quadratic 
instead of linear rise of the potential as a function of the distance. The last term, V4, is represented in Fig. 1(c) and corresponds 
to a connected flux tube. It is given by a Steiner tree, i.e., it is minimized by varying the location of the Steiner points si and S2. 
The choice of this potential is inspired by Refs. |[3j[T2l[l3][T4l, and has been discussed in the context of lattice QCD ||T5l[T6l . 




(a) Baryons (b) Tetraquarks: flip-flop (c) Tetraquarks: Steiner 



FIG. 1: Generalization of the linear quark-antiquark potential of mesons to baryons (left) and to tetraquarks, where the minimum is taken of 
the flip-flop (center) and Steiner tree (right) configurations. 
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The four-body problem in quantum mechanics is notoriously difficult. For instance. Wheeler proposed in 1945 the existence 
of a positronium molecule (e+ , e"*" , e~ , e~ ) which is stable in the limit where internal annihilation is neglected, i.e., lies below 
its threshold for dissociation into two positronium atoms. In 1946, Ore published a four-body calculation of this system lUTll 
and concluded that his investigation "counsels against the assumption that clusters of this (or even of higher) complexity can be 
formed". However, in 1947, Hylleraas and the same Ore published an elegant analytic proof that this molecule is stable [18J . It 
has been discovered recently [ 1^ . 

Similarly, the above model (j2]), in its linear version, was considered by Carlson and Pandharipande, who entitled their paper 
II20I "Absence of exotics in the flux tube model", i.e., did not find stable tetraquarks^ However, Vijande et al. |21| used a 
more systematic variational expansion of the wave function and in their numerical solution of the four-body problem found a 
stable tetraquark ground state. Moreover, unlike |20 |, they considered the possibility of unequal masses, and found that stability 
improves if the quarks are heavier (or lighter) than the antiquarks, in agreement with previous investigations (see, e.g., |21 1 for 
Refs.). 

It is thus desirable to check whether this minimal-path model supports or not bound states. The present attempt is based on an 
upper bound on the potential, which leads to an exactly solvable four-body Hamiltonian. 
With the Jacobi vector coordinates 

X = V2-vi , y = W4 - w 3 , z = , (3) 

and their conjugate momenta, the relative motion is described by the Hamiltonian 



H=% + ^ + f + U{x,y,z), (4) 
M m 4/i 



where /i, given by /i^^ m^^ + M^^, is the quark-antiquark reduced mass. Using the scaUng properties of H, one can set 
m — 1 without loss of generality. 

The simplest bound on the potential U is 

f/<^4< ||x|| + ||y|| + ||z|i , (5) 

as the tree with optimized Steiner points si and S2 is shorter than if the junctions are set at the middles of the quark separation 
V1V2 and antiquark separation V3V4. This leads to a separable upper bound for the Hamiltonian 

H<H'=^ + \\x\\+pl + \\y\\ + ^ + \\z\\. (6) 

Now, the ground state eo of + ||a;|| corresponds to the radial equation —u"{r) + ru{r) — eau{r) with u{0) = u{oo) = and 
is the negative of the first zero of the Airy function, eo = 2.3381 ... By scaling, the ground state of ap^ + P\\x\\, with a > 
and /? > is a^^^P'^^^eQ. Thus the lowest eigenvalue of H' is 



(7) 



with fj, = M/{1 + M). By comparison, the threshold of {QQqq) is made of two identical (Qq) mesons, each governed by the 
Hamiltonian h = /{2fi) + \\r\\, where p is conjugate to the quark-antiquark separation r. Thus the threshold energy is 

£;,h = 2eo(2^)-i/3^ (8) 

and it is easily seen that £" > E^h for any value of the quark-to-antiquark mass ratio M, i.e., the bound (jsj cannot demonstrate 
binding. 

A better bound will be proved below. If there is a genuine Steiner tree^ linking the quarks to the antiquarks, then 

^4<^(||x|| + ||y||) + ||z||. (9) 

But if V4 is not associated to a genuine Steiner tree, this inequality is often violated. Consider for instance a rectangular 
configuration with ||fiW2|| = 11^3^^411 ^ ll'^'i'^^sll — ll''^2'f^4|| (in this case the mathematical Steiner tree problem would require 



The authors used a relativistic form of kinetic energy and considered also the possibility of short-range corrections, but this seemingly does not affect their 
conclusion. 

^ This will be made more precise in the proof given in Appendix 
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a Steiner point linking vi and v^, another Steiner point linking V2 and 1)4, but the corresponding fluxes are not permitted by the 
color coupling in QCD), then ||z|| ^ and V4 ||a;|| + || so (j9]l does not hold. 
However, it will be shown that 



C/<^(lkll + ||2/||) + 11^11, 
for any configuration of the quarks and antiquarks, i.e., for any x, y and z. Then the ground state of H is bounded as 



(10) 



E<E" = eo 



1/3 



/3 



'1/3 



(11) 



As shown in Fig.[2j this bound E" significantly improves the previous one, E'. It is easily seen than E" becomes smaller than 
Ett for very large values of the mass ratio, more precisely for M > 6402, and thus that the tetraquark is bound at least in this 
range of AI. The numerical estimate of f2T] actually indicates stability for all values of M, even M — 1. 




M 
1000 



FIG. 2: Simple bound E' (Eq. (jv}, dotted line) and improved 
upper bound E" (Eq. 1 1 1 1, solid line) on the tetraquark ground- 
state energy as a function of quark-to-antiquark mass ratio M. 
Also shown is the threshold energy _Eth (Eq. [8j, dashed line). 
The energies are in units of eo, the ground state of — A + ||r||. 



To summarize, we obtained an analytic upper bound on the ground state energy of tetraquarks systems with two units of open 
flavor, {QQqq), using a model of linear confinement inspired by the strong-coupling regime of QCD. The key is an inequality 
on the length of a Steiner tree with four terminals. The bound confirms a recent numerical investigation, in which this potential 
was shown to bind these tetraquarks below the threshold for dissociation into two mesons. It remains to investigate whether this 
stability survives refinements in the dynamics, such as short range corrections, spin-dependent forces, etc. 

It is our intention to extend this investigation to the case of the pentaquark (one antiquark and four quarks) and hexaquark 
configurations (six quarks), which have been much debated in recent years. 
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for useful correspondence. 



Appendix: Results on the Steiner problem 



Before deriving Eq. ( 10 1, let us review some basic properties of the string potentials V3 and V4. 



Three terminals The three-point problem is very much documented in textbooks 1221 |23l l24l [25] |251 . Let be the 

triangle, with side lengths oi — \\v2V3\\, . . . and angles ai = ZV1V2V3, etc. The problem of finding a path of minimal length 
llsTJill + IISV2II + lls'^sll linking the three vertices has been solved by Fermat and Torricelli. See, e.g., [22|. The result is the 
following: if one of the angles, say ai, is larger than 120°, s coincides with vi, otherwise each side of the triangle is seen from 



s with an angle of 120°. The Steiner point s is thus at the intersection of three arcs of circles, see Fig. 3(a) 



The three-terminal problem is also linked to Napoleon's theorem, which states that if one draws external equilateral triangles 



on each side, V1V2W3, viv^wi and viViW2, the centers of these triangles form an equilateral triangle (dashed lines in Fig. 3(b) 1, 
a nice example of symmetry restoration. The junction s is just the intersection of viwi, V1W2 and v^. Note that \\svi \\ + \\sv2\\ — 
IISW3II, and similar relations, and thus the potential is simply V3 — |jwiWi|| — \\v2W2\\ = Wv^w^W- 

The point w-s and its symmetric with respect to V1V2, form the toroidal domain associated to the subset {vi,V2}- The length 
of the minimal Steiner tree is the maximal distance between V3 and the domain {ws, t^}. 
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From the above properties, one can estimate the string potential in a closed form. If ai > 120°, then V3 = 02 + 03, and 
similarly for large a2 or a^. Otherwise, V3 — Yl £i, where £i — ||stii||. From s = 120° in the triangle SV2V3, af = ^2 + ^3 + ^2^3, 
and by summation 

2{el +£l + £l) + {£i£2 + 4^3 + £3£i) = al + al + al. (12) 

Now, £i£2 being four times the area of the SV1V2 triangle, the second term in the above equation is four times the whole area of 
^^11^2^3, which is given by the Henon theorem. Altogether, in the case of a genuine Steiner tree ITJ 



V3 = ^1 + ^2 + 4 = V a? + ^2 + «3 + VH"-i + «2 + a3)(ai + ^2 - a3)(ai - 02 + a3)(-ai + 0.2 + 03) , (13) 
which can be computed quickly. 



The planar tetraquark problem. For the four-point problem, there are many special cases, which can be treated by inspec- 
tion. If, for instance the quark i>2 is on the back of vi, as in Fig. 4(a) the problem reduces to the Steiner problem for {wi, 113, V4}. 
Another special case is shown in Fig. 4(b) where the quarks are close to the antiquarks. For the standard Steiner problem of 
geometry, the solution would correspond to the Steiner tree shown as a dotted line, with a Steiner point S3 linked to vi and U3 
and another one, S4, linked to V2 and 114. This is not allowed by the different color properties of quarks and antiquarks, hence 
our best tree, shown as a solid line, has only one junction. But in estimating the potential U of Eq. (|2]l for this configuration, the 
minimum is the flip-flop term dis + ^24- 

Let us turn to the case of a genuine Steiner tree (wiW2)siS2(w3f4) as in Fig.js] The string of Fig. [T(c)| is minimized with respect 
to si and 82- Hence for fixed S2, it assumes the Fermat-Torricelli minimization of V1V2S1, a well-known iteration property of 
Steiner trees. Hence ZviSiV2 = 120° and V1V2 is the bissector of ZviSiV2 and passes through the point W12 which completes 
an equilateral triangle wiU2U'i2 in the quark sector. Similarly, it also passes through u'34 which makes V3V4W34 equilateral in the 
antiquark sector. 

The junction points si and S2 are just the other intersections of the straight line W12W34 with the circumcircles of 'i;it;2?«i2 
and V3V4W34, as shown in Fig.|5] There is a possible ambiguity about on which side si or S2 should be, but this is easily solved 
by the requirement that the total length of the string is minimum. Crucial is the observation that V — |l'ii'i2W34||, so that the 
determination of the Steiner points si and S2 is not required to compute V4. 

A variant is that is ti2 is the symmetric of W12 with respect to V1V2, the set {wi2ti2} is the toroidal domain associated to the 
quarks, and similarly {11134134} for the antiquarks, the length of the Steiner trees is the maximal distance between these two sets. 
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(a) Large angle invi (b) Elongated rectangle 



FIG. 4: Examples of special configurations. Left: one junction coincides with vi. Right: the two junctions merge (the dotted gray line 
corresponds to the Steiner tree if the four points Vi play the same role, unlike the tetraquark problem with quarks and antiquarks having 
conjugate colors. ) 




FIG. 5: Construction of the minimal string in the planar case. 

This construction, which is a special case of the Melzak's algorithm IZTl . leads to a very easy computation. If each vector Vi is 
identified with its affix (complex number) Vi, etc., then those of Wi2 and k;34 are easily deduced, for instance Wi2 = —j'^vi ~jv2 
or —jvi — j'^V2 (depending on which side is W12), if one uses the familiar root of unity j — exp(2i7r/3). Once W12 and k;34 
are determined, V — ||wi2W34||. If one wishes to locate the Steiner points, it is sufficient to remark that wi2S2.u'i2'i«34 — 
\\wi2C34\\^ — and W34S1.W34W12 — ||w34Ci2||^ — 7-12, where C12 is the center of the circle z;i-(;2Wi2 and ri2 = rfi2\/3/2 its 
radius and C34 and are defined similarly in the antiquark sector. 

The spatial tetraquark problem In general, the four constituents do not belong to the same plane. The minimum is 
achieved for V1V2S1S2 coplanar, and V3V4S1S2 also coplanar, but in a different plane. The toroidal domain to which the point 
W12 belongs is the Melzak circle, of axis V1V2 and radius ri2 = ||i;it;2||\/3/2, and similarly for 1034 in the antiquark sector. The 
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straight line W12W34 has to intersect these two circles as well as the lines V1V2 and V3V4. The problem consists of constructing 
such a straight line. 

The reasoning can be made on Fig. [5] if one imagines that U3W4S2 is not coplanar to V1V2S1. As stressed in ||28l . the key is to 
determine p and q, the intersections of S1S2 with V1V2 and V3V4, respectively. In this paper, the following coupled equations are 
derived 

niy/h'^ + x'l SIT? (j) + ri2V cos (j) + xl sit? (j) + r34V cos (p 

- - ^ - - ^ , (14) 



+ x"^ SIT? (j) ri2 + + x"^ SIT? 4> 



for the abscissa Xp of p along V1V2 and Xq of q along v^v^. These abscissas are from the common perpendicular uv to V1V2 
and vy,Vi (u £ t;iz;2and v E ^3^4), with \\uv\\ — h, \\uh\\ — m and ||wfc|| = n, where h is the middle of V1V2 and k that 
of V3V4. The equations ( [T4| i can be solved by iterations, with remarkably fast convergence. Once Xp and Xq, i.e., p and q, 
are determined, the Steiner points are determined by imposing they are on the circles witi2U'i2 and V3V4W34, respectively. For 
instance, if si — p + t{q — p), t obeys a second order equation-'. 

If one is interested only in the length of the Steiner tree and not in the position of the Steiner points, an alternative formalism 
consists of locating p through p ^ h + x {v2 ~ h) and q = k + y {v4 ~ k) . With this notation, the length of the tree is simply 



V4 — min 

x,y 



(15) 



which is easily minimized by varying x and y. The minimisation is equivalent to solving the coupled equations 



X = ^/3 + x'^ ^^^n'!i'^ II ' y = V^ + y'^ u '"^^tT II ' 

which expresses that W12, p, si, S2, q and W34 are collinear. These equations are easily solved by iteration or any other means. 
We believe that, besides checking the particular cases with large angles or a single Steiner point, the fastest computation of 



the connected four-quark potential consists of minimising (15 1 or solving ( 16 1. We expect a dramatic improvement in computing 
time from the above algorithm. 

However, it is aesthetically appealing to attempt a further reduction of the number of variables to be determined numerically, 
and to provide an almost analytic estimate of the interaction as a function of the coordinates of the quarks and antiquarks. 
Finding V4 = ||t«i2W34||, the maximal distance between the Melzak circles C12 and C34, is very similar to the problem of the 
minimal distance beween two circles in space, as addressed e.g., in |29;,,30J. Neff I.29J has shown that with the help of Lagrange 
multipliers and Grobner type of elimination performed by computer-algebra sofware, the squared stationary distance V4 obeys 
an eighth-order polynomial equation whose coefficients are rational functions of the coordinates of vi, V2, V3 and V4. 

Eberly |30| showed that if m is associated to an angle 9 along C12, and n to cj) along C34, then imposing ||mn||^ to be stationary, 
results in two equations of the type 

aiCOsO + PiSine + 'f, = , i = l,2, (17) 



where a^. Pi and 7^ contain constants and terms linear in cos </> and sin Solving ( 17 1 as two linear equations, as if cos 9 and 
sin 6 were independent, and then imposing cos 9^ + sin 9^ = 1 gives an equation for cos and sin 0, which is transformed into 
an 8*^ order equation in cos cj) 



It is slightly faster to rewrite ( 17 1 using t — tan(0/2) and u — tan(0/2) as 

5,t^ + T],t + e, ^ , j = l,2, (18) 
where the coefficients are quadratic in u. The compatibility of two such equations is simply 

WiS,r,)W{7j,e)^W{6,ef , W{x,v) ^ xiy2 ~ X2yi , (19) 
and is directly a polynomial in u, of order 8. 



' There is a misprint in 1 28 1 whicli propagated in tlie numerical calculation given as an example. 
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FIG. 6: The confining potential V4 for tlie tetraquark system {V1V2VSV4) is the 
minimal length of the tree || + ||«2Si || + ||siS2 1| + ||s2ti3 1| + ||s2U4|| when si 
and S2 are varied. It is also the maximal distances between the circles C12 and C34, 
i.e., the distance «;i2iA'34. The Melzak circle C12 is centered at the middle of vi and 
V2, has 111112 as axis and a radius \\v1V2 
the antiquark sector. 



|\/3/2, and C34 has analogous properties in 



Proof of the inequality ( [TO) If we have a positively oriented edge from si to S2, i.e., the Steiner tree is non degenerate, then 
we have 

V4<{\\x\\ + \\y\\)V3/2+\\z\\=B 

using Melzak circles. 

However the bound required is for U — min {di^ + ^24, ^14 + ^23, V4}. So we want to confirm that 

U <{\\x\\ + \\y\\)V3/2+\\z\\^B 

is valid, regardless of whether V4 is a degenerate or non degenerate Steiner tree. 

We follow the variational method introduced in |31 1. The problem is formulated as a global optimisation problem as follows; 

Define L as the length of the formal Steiner tree spanned by the four vertices. This length is obtained from the distance between 
the farthest points on the two Melzak circles. In terms of the usual Steiner tree components, L = ||wisi|| + ||w2Si || ± |lsiS2|| + 
l!''^3S2|| + ||i'4S2||)- We get the positive sign for ||siS2|| if there is a real Steiner tree. On the other hand, if the Steiner vertices 
have interchanged position, so that on the line between the two farthest Melzak points, S2 is closer to the Melzak point for vi , V2 
than si, then we have the negative sign for ||siS2||. So we can construct a formal tree on the six vertices wi, W2, f3, va, si, S2 
where the edge joining the two Steiner vertices is 'negatively oriented'. 

Now it is easy to see that L < {\\x\\ + ||y||)\/3/2 + ||z||. Soify = L then the desired inequality follows trivially. So we only 
need to consider the situation where L < V, i.e the Steiner tree is formal rather than a real Steiner tree. Now by the inequality 
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above, if either of dia + (i24, du + ^23 is not larger than L, then clearly the required inequality follows. So we only need to 
consider the case when di3 + d24 > L and du + (^23 > L. 

We can parametrise the points Vi,V2, wa, W4 by the numbers ||tiiSi||, ||u2Si||, ±||siS2||, ||i'3S2||i 1^45211- (It is easy to see that 
these four points are determined up to rotation, translation by five parameters.) By rescaling, we can assume that the sum of these 
five numbers is 1, without loss of generality for the inequality. It is easy to see that all the numbers are then bounded so the domain 
becomes compact. So we seek a maximum of the ratio of R = mm{di3 + 0^24! c?i4 + ^23} and (||a;|| + ||j/||)\/3/2 + ||z|| = B 
over this domain. 

Now suppose that we rotate the triangles V1V2V3 and V1V2V4 around an axis line through viV2- Clearly we can think of one 
triangle as being fixed and the other as moving relative to the first one. The quantity R does not change by this rotation, but 
obviously B does. Hence a maximum of the ratio R/B corresponds to a minimum for B under such a rotation. 

Now an elementary argument shows that such a minimum for B occurs for the configuration being planar, i.e when the vertex 
V4 moves into the plane of vi,V2,V3. Now assume that some initial configuration satisfies R/B > 1 and the Steiner tree is 
formal rather than real. As the triangle V1V2V4 rotates around an axis line through V1V2, it is easy to see that the two Melzak 
circles move apart. At some intermediate point, if they cross, then we find that the Steiner tree changes from being formal to 
being real. At this intermediate point, it is trivial to see that R/B < 1. But this is impossible, since we have initially R/B > 1 
and R/B is increasing, since B is decreasing and R is fixed. 

On the other hand, if the Melzak circles never intersect, then this must be true for the planar configuration. So we would have 
such a configuration for which the Steiner tree is still formal but R/B > 1. It is elementary to prove that this is impossible. So 
this completes the argument. 
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